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Recent observations seem to suggest that our Universe is accelerating implying that it is dominated 
by a fluid whose equation of state is negative. Quintessence is a possible explanation. In particular, 
the concept of tracking solutions permits to adress the fine-tuning and coincidence problems. We 
study this proposal in the simplest case of an inverse power potential and investigate its robustness 
to corrections. We show that quintessence is not affected by the one-loop quantum corrections. 
In the supersymmetric case where the quintessential potential is motivated by non-perturbative 
effects in gauge theories, we consider the curvature effects and the Kahler corrections. We find that 
the curvature effects are negligible while the Kahler corrections modify the early evolution of the 
quintessence field. Finally we study the supergravity corrections and show that they must be taken 
into account as Q « mpi at small red-shifts. We discuss simple supergravity models exhibiting the 
quintessential behaviour. In particular, we propose a model where the scalar potential is given by 
V{Q) = A qo e?® . We argue that the fine-tuning problem can be overcome if a > 11. This model 
leads to ujq « —0.82 for f2 m w 0.3 which is in good agreement with the presently available data. 

PACS numbers: 95.35+d, 98.80.Cq 



I. INTRODUCTION 

Several observations seem to suggest that our present 
Universe is dominated by a type of matter with nega- 
tive equation of state, ujq = Pq/pq < 0. The first type 
of observations leading to this conclusion is the recent 
measurements of the relation luminous distance versus 
redshift using type la supernovae jlj . The interpretation 
of the data are usually made under the assumption that 
the unknown fluid is a "true" cosmological constant A. 
Unfortunately, the results are degenerate in the f2 m — 
plane and it is difficult to draw a conclusion on the basis 
of these measurements only. The situation changes dras- 
tically if one includes in the analysis a second type of ob- 
servations: the measurements of the CMB anisotropics. 
In this case the degeneracy can be removed ^ and one is 
led to the conclusion that the matter with negative equa- 
tion of state would contribute by 70% to the total energy 
of the Universe, the remaining 30% being essentially Cold 
Dark Matter ensuring that the Universe is spatially flat, 
ilo = 1 , in agreement with the standard inflationary sce- 
nario. This conclusion can be traced back to the fact that 
many CMB experiments show a high amplitude of the 
first Doppler peak located at £ « 260. For example, this 
is the case for the experiments Saskatoon fl, PythonV 
j§ or TOC097 f§. The addition of a fluid with negative 
equation of state has for consequence that the Integrated 
Sachs- Wolfe effect reinforces the scales I w 200 — 300 
and increases the peak to values compatible with the er- 



ror bars of these experiments. Moreover, the position 
of the peak informs on the value of f^o where fio is the 
ratio of the total energy density to the critical energy 
density. In the simplest case the position is predicted to 
be ^Doppler ~ 220/\/r2o. In the case where a fluid with 
negative equation of state is added the peak is shifted 
towards bigger values of t. This seems to be the case for 
the experiments cited above. 

Recently, another method was proposed in order to re- 
move the degeneracy between f2 m and fl\ using large- 
scale peculiar velocities data ||. These data provide 
constraints mainly on J7 m and are almost independent 
of f^A- Combined with the measurement of the relation 
luminous distance versus redshift, they select a region in 
the Q m — 57a plan which is compatible with the results of 
Ref. J2|. It is remarkable that, although of different na- 
ture, these experimental data converge towards the same 
conclusion. 

This raises the issue of the physical origin of this fluid 
with negative equation of state. A useful indicator of 
the physical nature of this fluid is the value of ujq . Re- 
cent constraints Jml indicate that —1 < ujq < —0.6 
whereas in Refs. g|| a value such that — 1 < ujq < — 0.8 
is favoured. The case ujq = —1 corresponds to the 
existence of a "true" non-zero cosmological constant. 
This cosmological constant has then to be explained 
by current particle physics scenarios. In particular one 
has to face the task of explaining an energy scale of 
w 5.7ft 2 x 10~ 47 GeV 4 , i.e. a value far from the natu- 
ral scales of particle physics. Therefore, although per- 
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fectly compatible with the presently available data, this 
hypothesis runs into theoretical problems since it seems 
easier to explain a vanishing cosmological constant (by 
some yet unknown fundamental mechanism maybe com- 
ing from quantum gravity or string theory, see Ref. |10|) 
than finding a reason for a tiny (in comparison with 
the high energy physics scales) contribution. In a cer- 
tain sense the measurements described above render the 
"quantum" cosmological constant problem worse than 
before. 

Recently, another explanation, named quintessence, 
has been put forward in Refs. |fl]]. Quintessence is an 
alternative scenario with a homogeneous scalar field Q 
whose equation of state is such that — 1 < ujq < 0. In this 
scenario, the missing energy density is due to this scalar 
field. Let us note that this explanation allows to come 
back to the situation where there is a vanishing cosmo- 
logical constant. However the quintessence scenario does 
not solve the "quantum" cosmological constant problem. 

Quintessence has to adress a certain number of ques- 
tions. First of all one must make sure that the fine tuning 
problem of the cosmological constant does not reappear 
in a different guise. One must also solve the "coinci- 
dence problem" , i.e. understand why the quintessential 
field begins to dominate now. Another conundrum is 
to try to justify the presence of such a field from the 
Particle Physics point of view. The answers to these 
questions strongly depend on the form of the potential 
V(Q). For example, if one chooses a potential of the 
form V(Q) = (l/2)m 2 Q 2 then one cannot avoid to fine 
tunc the value of the mass to an extremely small num- 
ber jL2| . The problem is then similar to the case of the 
cosmological constant. 

However, the problems described previously can be 
adressed if one considers the following potential |l3|] : 



V(Q) = 



A 



a+4 



(1) 



where a > and A are free parameters. This poten- 
tial possesses remarkable properties. The equations of 
motion have an attractor solution called in Ref. G3| the 
"tracking field". The initial conditions can vary by 100 
orders of magnitude leading to the attractor in all the 
cases. Since the present value of Q ss mp\ on the attrac- 
tor, one has that A w (f^pcm^) 1 ^ 4 ^ « 10 6 GeV for 
a = 6 where p c = 8.1h 2 x 10 _47 GeV is the present value 
of the critical energy density. This value is not in con- 
tradiction with usual high energy scales. Moreover, one 
can hope to justify the form of the potential given in Eq. 
(|l|) from high energy physics p^-p^]. Finally it should 
be noted that, in principle, it is possible to distinguish 
quintessence from a cosmological constant since one has 
in general ujq ^ — 1. 

The aim of this paper is to study the robustness of the 
concept of tracker solutions. In section II we quickly re- 
view the main properties of the tracking solutions. Then 
in section III, we analyze whether the nice properties of 



the tracking field are affected by the quantum corrections 
to the potential given by Eq. ([l]) at the one loop level in 
the case where the underlying model is not supersymmet- 
ric. We show that the quintessential scenario is robust 
against these corrections. In section IV, we turn to the 
study of the SUSY models. We argue, as already noted in 
Ref. , that potentials given by Eq. ([!]) naturally arise 
in the context of supersymmetric gauge theories where 
certain flat directions are lifted by non-perturbative ef- 
fects. We study the phenomenology of these models for 
which the quantum corrections to the superpotential au- 
tomatically cancel out. As quintessence requires a value 
of the second derivative of the potential of the order of 
the scalar curvature of the universe, we also study the 
corrections due to the fact that the fields live in a curved 
spacetime. The curvature effects are evaluated at the 
one-loop level and shown to preserve the tracker field 
properties. Finally in the supersymmetric case one can 
take into account the effect of the corrections to the ki- 
netic terms of the quintessence field. In particular in the 
low energy description of the supersymmetric gauge the- 
ory the Kahler potential receives corrections suppressed 
by the gaugino condensation scale. We show that this 
leads to difficulties for the supersymmetric models of the 
tracker potential. In section V, we turn to the study of 
SUGRA models of quintessence. We emphasize that such 
models are the most physical ones since at the end of the 
evolution the field is on tracks which implies that its value 
today is Q w mp\. We analyse the SUGRA corrections 
to the inverse power law potential and show that they 
lead to inconsistencies due to the possible negative val- 
ues of the potential. To remedy this situation we propose 
a supergravity scenario where the potential is guaranteed 
to remain positive. We apply this framework to the case 
of the heterotic string where the role of the quintessence 
field is played by the string moduli. Indeed the moduli 
are famous for leading to run-away potentials as expected 
for quintessence. We find that the resulting potential is 
exponentially decreasing, a case already studied in the 
literature which fails to give the appropriate energy den- 
sity. We eventually present a toy model where the inverse 
power law results from the SUGRA potential. We end 
with the conclusions presented in section VI. 



II. TRACKING SOLUTIONS 

In this section we quickly review the main properties 
of the tracking solutions as explained in Refs. p3[ . 

The Universe is described by a spatially flat 
Friedmann-Lemaitre-Robertson- Walker (FLRW) space- 
time whose metric can be written as: ds 2 = — dt 2 + 
a 2 {t)dx 2 . 

We assume that the matter content of the Universe 
is composed of five different fluids: baryons, cold dark 
matter, photons, neutrinos and the quintessential field 
Q. The energy density of Baryons and cold dark matter 
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evolves as p m = p c f2 m (l + z) 3 where z is the redshift. 
The equation of state is p m — which is equivalent to 



0. Observations indicate that fL 



-+- Q c 



0.3. Photons and neutrinos have an energy density given 
by Pr = Pc^r(l + z) A . The equation of state is given 
by u> r = 1/3. The contribution of radiation is negligible 
today since = n y + Cl v ps 1CT 4 . Finally, the fifth 
component is the scalar field Q. Its equation of state 
is characterized by uj q = [\Q 2 - V{Q)]/[\Q 2 + V{Q)\ 
where a dot represents a derivative with respect to the 
cosmic time. A priori, ujq is not a constant and is such 
that —1 <ujq < 1. Since the Universe is supposed to be 
spatially flat, we always have f2 m + f2 r + fig = 1 which 
leads to J7q « 0.7. In the following, we will denote the 
dominant component in the energy density by pb so that 
during the radiation dominated era we have pb = p r and 
during the matter dominated era, pb — p m - A similar 
notation will be used for wb- 

The evolution of the scale factor is governed by the 
Fricdmann equation: 



H 2 



:(pr, 



pq), 



(2) 



where k = 8irG/c 4 = 87r/mpj in the Planck system of 
units. The evolution of the scalar field is given by the 
Klein Gordon equation: 



Q + 3HQ + V'(Q) = 0, 



(3) 



where a prime denotes the derivative with respect to Q. 

The inverse power law potential was first studied in 
Rcf. If one requires that, during the radiation dom- 
inated era, the energy density of the scalar field be sub- 
dominant (this is necessary for not being in conflict with 
the Big Bang Nucleosynthesis), i.e. pq <C pb, and red- 
shift as pq oc a~ 4a /( a+2 ) then one is automatically led to 
the potential of Eq. (|l|). This was the original motiva- 
tion of Ref. |51| for considering the potential ([!]). In that 
case, it is possible to find an exact solution to the Klein 
Gordon equation for which Q oc a 4 ^ a+2 \ One can show 
that this solution is an attractor . Then, if one follows 
the behaviour of the scalar field during the matter dom- 
inated era (i.e. for pb = p m ) with the same potential, 
one can show pl| that Q oc a 3 /( Q+2 ' is an exact solu- 
tion which is still an attractor. For this solution, one has 
Pq oc ct~ 3Q /( a + 2 ). The previous results are equivalent to 
saying that the attractor is given by: 



d 2 V(Q) 9a + l 



(l~u 2 )H 2 , 



(4) 



dQ 2 2 a 
during both the radiation and matter dominated epochs^ 



* Establishing this relation from Ref. |l8| the factor 1/2 is 
no longer present because the definition of the potential in 
that paper differs by this factor from the definition adopted 
in Ref. [JL8| and in the present article 



We can re-write the parameter ujq as ujq — (awe — 
2)/(a + 2). Since pq redshifts slower than pe, the scalar 
field contribution becomes dominant at some stage of the 
evolution. 

As shown in Ref. pa], this scenario possesses impor- 
tant advantages. Firstly, as already stated in the in- 
troduction, one can hope to avoid any fine-tuning. In- 
deed if the scalar field is on tracks today and begins 
to dominate and if, in addition, we require Qq « 0.7 
then A w 4.8 x 10 6 GeV (for a = 6), a very reason- 
able scale from the High Energy Physics point of view. 
Secondly, the solution will be on tracks today for a 
huge range of initial conditions. If one fixes the ini- 
tial conditions at the end of inflation, z = 10 28 , the al- 
lowed initial values for the energy density are such that 
10~ 37 GeV 4 < p Q < 10 61 GeV 4 where 10" 37 GeV 4 is ap- 
proximatively the background energy density at equal- 
ity whereas 10 61 GeV 4 represents the background energy 
density at the initial redshift. If the scalar field starts 
at rest, this means that 10~ 18 mpi < Qi < 10~ 2 mpi ini- 
tially. Thirdly, the value of ujq is automatically such that 
— 1 < uq < today. The precise value of ujq depends on 
the functional form of V(Q) and on the value of f2 m . 

In the following figures, we illustrate these properties 
for a = 6, ft m = 0.3 and Qi(z = 10 28 ) « 3 x 10- 18 m P1 
which roughly corresponds to equipartion at that red- 
shift, that is to say Oq^ ps 10~ 4 . Equations (^) and (||) 
are integrated numerically. The first figure represents the 
evolution of the energy densities throughout the radiation 
and matter dominated epochs. 
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FIG. 1. Energy density of radiation, matter and 
quintessence in function of the redshift starting from equipar- 
tition. Dotted line represents radiation, dashed line is matter 
and the full line is quintessence. 

The interpretation of these curves has already been 
given in Ref. |^3[ . The case presented here corresponds 
to an "overshoot" according to the terminology of that 
reference. First the scalar field rolls down the potential 
such that its kinetic energy dominates and Q oc a~ 3 . 
Then, the field freezes to some value Qf. And finally, it 
joins the attractor. 
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In the next figure, the evolution of the equation of state 
for the same model is displayed. 
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FIG. 4. Equation of state for quintessense starting at 
'; » 0.2 x 10" 4 . 



FIG. 2. Equation of state in function of the redshift start- 
ing at equipartition. 

The value of ojq today for this model is found to be 
ujq « —0.4. Therefore it is clear that this case cannot be 
considered as a realistic case but rather as a toy model. 

In order to illustrate the insensitivity to the initial con- 
ditions, the two following figures show the same case as 
previously but with an initial value of the scalar field 
given by Qi w 0.2 x 10~ 4 m P i. 
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FIG. 3. Energy density of radiation, matter and 
quintessence in function of the redshift starting at 
Q; « 0.2 x 10~ 4 . 

This case corresponds to an "undershoot" . One can see 
that the field starts directly from its frozen value. Finally, 
the evolution of the equation of state is displayed. It is 
apparent that it leads to the same cosmology today with 
ujq « -0.4 



In the next section we study the influence of the quan- 
tum corrections to the potential (|l|) on the properties 
described in this section. 



III. QUANTUM CORRECTIONS TO NON SUSY 
MODELS 

At the classical level we have chosen a potential given 
by Eq. ([!]). However, it is a generic effect that this po- 
tential will be modified when quantum corrections are 
taken into account. In this section, we only study the 
one loop corrections. These types of corrections auto- 
matically cancel out when the model is supersymmetric. 
Other corrections such as the corrections due to curva- 
ture effects and to the kinetic terms will be studied in 
the next sections. 

The modified potential reads U$ BjJ: 



V(Q) 



Q a 4 


A4 l 


A 2 


TO 4 


( m 2 
l ln A2 " 


r 


+ 32. 2 





2A2 



m 2 A 
32tt 2 



(5) 



where A is an effective cut-off already defined and fi is 
the natural energy scale of the theory. This expansion is 
obtained by calculating the one loop Feynman diagrams. 
This amounts to evaluating the integral J d 4 pln(p 2 +m 2 ) 
properly regularized by A. This choice has been made be- 
cause A turns out to be the natural cut-off in the physical 
models considered in this paper Q. The energy scale 
/i appears in the renormalization conditions. It turns 
out that its precise value is not important for our pur- 
pose since it only appears in the logarithm. The effective 
"mass" m 2 is equal by definition to: 



2 <?V(Q) A 4 
m = — — = <x(a + 1) 



dQ 2 



Q 



a+2 ' 



(6) 



The second term in Eq. (|5|) does not depend on the field 
Q. This term will contribute as a cosmological constant. 
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Of course all the other fields in the Universe also give 
contributions to the cosmological constant. It is hoped 
that, by some unknown mechanism, the total contribu- 
tion vanishes, see the introduction. This assumption is in 
the spirit of the quintessential models in which there is no 
need of a cosmological constant in the Einstein equations. 
For all these reasons we will not consider the second term 
in Eq. (|j) in what follows. 

Introducing the expression giving m 2 into the formula 
of the corrected potential, one finds: 



V(Q) 



A 4+Q 



t(a + l)A 6+Q 1 



Q a 

q 2 (q + l) 2 A 8 + 2 ° 

2q+4 



32tt 2 
1 



32tt 2 



Q 



In 



Q a+2 
a(a + 



1)A 2 



)a+2 



(7) 



We see that the functional form of the potential is no 
longer the same. 

We now need to estimate the orders of magnitude of 
the corrections to sec whether they can be important. 
As an example, let us consider the case a = 6 for which 
A 4.8 x 10 6 GeV. The first change is that, now, we must 
have Qi > 4 x 10~ 15 mpi initially in order that pq < pb 
at z = 10 28 . Interestingly enough this constraint comes 
from the last term in Eq. (Q) which is dominant at this 
redshift. This means that there exists a region for which 
the quantum corrections are more important than the 
unperturbed potential. As expected, this happens in the 
early Universe, at high energy. A quick estimate en- 
ables to show that quantum corrections are dominant if 
Qi < 2 x 1(T 13 topi, i.e. for 10 25 GeV 4 < p Q < 10 61 GeV 4 
initially. Therefore, among the 95 orders of magnitude in 
which the initial energy density of quintessence can vary, 
35 of them are dominated by quantum corrections includ- 
ing the most physical case of the equipartition for which 
we have Qi as 6x 10~ 15 topi corresponding to Q,Qi ss 10~ 4 . 
We conclude that, a priori, quantum corrections must be 
taken into account in any realistic model of quintessence. 

However, we find numerically that the final value of pq 
and loq is the same with and without quantum correc- 
tions even if we start from equipartition. As this conclu- 
sion is not changed if one considers other values for a, we 
have demonstrated that quintessence possesses another 
remarkable property: it is stable against quantum cor- 
rections. In fact the evolution of both the energy density 
and the equation of state with and without the quantum 
corrections is the same during all the cosmic evolution. 
This is due to the fact that at the beginning of the evo- 
lution the field rolls down the potential very quickly and 
leaves the region where quantum corrections are impor- 
tant in a very short time. 

In conclusion we have shown that non SUSY models of 
quintessence are stable against one loop quantum correc- 
tions to the effective potential. This property is generic 
and does not depend on the precise value of a. 

In the next section, we start examining SUSY models 
of quintessence. In that case the quantum corrections 
studied in this section automatically cancel out. 



IV. SUSY MODELS 

As seen in the previous sections the quintessence field 
varies over a large range of values, as high as the Planck 
scale. It is therefore compulsory to treat the quintessence 
behaviour within the framework of models encapsulat- 
ing the expected behaviour of high energy physics. We 
will use supersymmctric models. One of the advantages 
of these models is their stability with respect to quan- 
tum corrections. In particular it is known that the non- 
renormalisation theorem preserves superpotentials, there 
is only a wave function renormalisation. 



A. SUSY gauge theories 

The potential in 1/Q a which leads to quintessence has 
a natural supersymmctric origin. This was first shown 
in Ref. [Q. The aim of this section is to generalise the 
models investigated in Ref. |l4| . This type of potential is 
generated at low energy E < A in supersymmetric gauge 
theories due to non-perturbative effects along flat direc- 
tions of the scalar potential. In this section we present 
the general features of supersymmetric gauge theories. 
More details can be found in [ P4j23f| 

Before dealing with the general case let us recall the 
most famous example constructed with the gauge group 
SU(N C ) and Nf quarks and antiquarks (z^,z^), i — 
I . . . Nf and a — 1 . . . N c , in the fundamental and an- 
tifundamental representations of the gauge group. The 
dynamics of this model is governed by the renormali- 
sation evolution of the gauge coupling constant. For 
b = 3N C — Nf > the gauge coupling constant becomes 
strong at low energy. In this infrared regime it is rele- 
vant to study the configurations of the scalar components 
of the quarks and antiquarks which have zero energy. 
The potential Vr> = g 2 D 2 /2 vanishes when the Z)-terms 
D a = z^T^zjj - z^T^z^p vanish leading to: 



6 a 



(8) 



The manifold of solutions of these equations is called the 
moduli space of the gauge theory. The moduli space is 
in one to one correspondence with the gauge invariant 
polynomials I a via the equation 



Qja 



(9) 



In the case Nf < N c there is just one gauge invariant 
I a = M y = z z a z^ called the meson field. The low energy 
dynamics is expressed in terms of the meson field. As 
the gauge group becomes strongly interacting at low en- 
ergy E < A where A is the strong interaction scale, non- 
perturbative effects due to the condensation of gauginos 
lead to a non-zero superpotential for the meson field, 



W oc 



A 6/(iV c -iV f ) 



(det M) 



1/(N C -N { ) ■ 



(10) 
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This potential is deduced for Nf = N c — 1 by an instanton 
calculation and for all Nf using the decoupling technique. 
Let us now focus on the amplitude mode M 1 -? = Q 2 8 l i 
obtained from z l — Qz$, z 1 = Qz~q. Note that the case of 
two field directions has recently been investigated in Ref. 
jLjJ . Starting from a flat Kahler potential K = tr(zz^ + 
zz$) the classical Kahler potential at low energy becomes 
K = QQ* after normalizing tr (zqZq + zqZq) = 1. This 
yields the low energy scalar potential for the amplitude 
mode 



V(Q) = 



A 26/(iV -JV f ) 
Q2(Nc+N { )/(N c -Ns) ■ 



(11) 



From the previous equation, it is clear that the so far 
arbitrary coefficient a is now given by: 



a 



) Nc±N L 
J N C -N f ' 



(12) 



which is always greater than two. This result had already 
been obtained in Ref. jL4|. The analysis carried out in 
this particular model is in fact typical of most supersym- 
metric gauge theories. We now show how this approach 
can be generalized. 

Consider a gauge group G and the matter fields z l in 
the representation i?j of the gauge group. Define the 
Dynkin index [i\ of a representation R\ to be 



tr(T a T fc ) = ££J' 



(i b 



(13) 



where the T a 's are the Hermitean generators of the rep- 
resentation Ri and 8 is the long root of the Lie algebra of 
G. For instance for G = SU(N C ), the previous results are 
retrieved using /if = 1 for the fundamental representation 
and /i a( jj = 2N C for the adjoint representation. The beta 
function determining the evolution of the gauge coupling 
constant depends on 2b = 3/i a dj — (J, where fj, = J2i i s 
the sum of all the Dynkin indices of the matter fields. The 
gauge theory is asymptotically free and strongly coupled 
in the infrared when b > 0. At low energy along the flat 
directions determined by the vanishing of the D terms 
the dynamics of the gauge theory is encoded in the prop- 
erties of the polynomial gauge invariants I a . When the 
Dynkin indices satisfy /i < /j, a dj the ring of gauge invari- 
ants is free. In terms of these gauge invariants the strong 
non-perturbative dynamics of the gauge theory generates 
a superpotential W(I) whose form is dictated by symme- 
try arguments. This superpotential when expressed in 
terms of the original matter fields z % reads 



W{Q) oc 



A 6 /(Madj-M) 



2/(/iadj-M) 



(14) 



As for the SU{N C ) case one is interested in the amplitude 
mode Q whose classical Kahler potential is flat. This 
leads to the scalar potential: 



V(Q) = 



^26/(Madj-M) 
Q2(Ai adj +M)/(Madj-M) ' 



(15) 



where the exponent a is now given by the following ex- 
pression: 



, Madj 
Madj ^ At 



(16) 



which is always greater than two. 

Therefore, the inverse power potential is a generic pre- 
diction of supersymmetric gauge theories. Nevertheless 
it is dependent on the hypothesis that the Kahler poten- 
tial is flat. This is plausible when Q > A but needs to be 
reconsidered below this scale. This leads to the danger- 
ous Kahler corrections which will be studied later. It has 
also been assumed that spacetime was flat. Curvature 
corrections are studied in the next section. 



B. Curvature corrections 

In the previous section we have treated the globally 
supersymmetric case as if the spacetime was not curved. 
These curvature effects are usually neglected as the typi- 
cal scale H is too small compared to the particle physics 
scales (recall that Hq ps 2.1h x lCP 42 GeV). Concern- 
ing quintessence this assumption has to be carefully 
checked as Eq. (|j) implies that the effective mass of 
the quintessence field is also of order H . This requires to 
study in a painstaking fashion the effects of curvature on 
global supersymmetry. A general argument is presented 
in an appendix while an explicit calculation is given here. 

We assume that the quintessence field belongs to a 
chiral supermultiplet, i.e. a complex scalar field and a 
Weyl spinor. We shall first examine the case of a free field 
showing that global supersymmetry is broken explicitly 
by curvature effects in a FLRW spacetime. This is due to 
the fact that the bosonic particles are created from the 
gravitational background whereas this is not the case for 
the fermions since they are conformally invariant. In the 
interacting case the breaking of supersymmetry by the 
curvature follows from a supergravity argument. In that 
case we evaluate the effective action at the one-loop level 
using the zeta function technique. 



1. Complex scalar field 

The action of the complex scalar field (pix^) is given 
by the following expression: 



S v = - ^xj—g9^d^d v v\ 



(17) 



It is convenient to separate the real and imaginary parts 
of the field and to write: 



1 



<P = +^2), 



(18) 



G 



where (fi, i = 1,2, are now real scalar fields. Since the 
spatial sections are flat, the fields can be Fourier decom- 
posed. It is convenient to perform this decomposition 
according to: 



^,x) = ^^/d»^ l k)e** 



(19) 



where we have extracted a factor 1/0(77) m the- time de- 
pendent amplitude of the Fourier component for future 
convenience. The Fourier components fJ,i(r], k) are such 
that (J,i(rj, k)* = Hi(r], — k) because the fields ipi are real. 
It can be easily seen that jUj(?7,k) obeys the equation of 
a parametric oscillator: 



k M»(»7) k ) 

a 



0. 



(20) 



This equation reduces to the equation of an ordinary har- 
monic oscillator in a Minkowski or radiation dominated 
Universe where a" = 0. 

When quantization is carried out the fields (pi become 
operators. In the canonical approach, the complex scalar 
field can be expressed as: 



0(i7, x) 



1 



1 



a(r?) (2tt) 3 / 2 
d 3 k 



4k0?) e 



— tk-s 



(21) 



where the operators djk are related to the operators of 
creation and annihilation Cik(rj) and cj k (?7) associated to 
the field operators ipi and satisfying the commutation 
relation [cjk(r?), ct k , (77)] = %<5(k — k') through the ex- 
pressions: 



Clk + «C 2 k 



V2 



it = 

- t 2k — 



-lk 



- IC. 



2k 



V2 



(22) 



The expression (^l|) can be used in order to express 
the Hamiltonian operator for the complex scalar field in 
terms of creation and annihilation operators. The result 
reads: 

H = \Y.J d3k ( fc (c*k4 + 4„ kCj ,_ k ) 

2 — 1 ^ 



■ a I t t \ 

1— (CikCi.-k - c! k c| _ k ) 



(23) 



The first term of the Hamiltonian represents the Hamil- 
tonian of a collection of harmonic oscillators. The sec- 
ond term represents the interaction between the classical 
background and the quantum field. It is proportional to 
the first derivative of the scale factor and therefore van- 
ishes in the Minkowski case. This term is responsible for 
the creation of (pair of) particles. If we start from the 
vacuum state |0) (no particle), then due to the presence 
of this interaction term, the state will evolve in a vacuum 
squeezed state |24]. 



2. Weyl spinor field 

The action of the Weyl spinor field is given by the 
following expression: 

1 



(24) 



Sip = - J / d 4 a:y^( M^T^^a 



In this expression E M are the Pauli matrices in a FLRW 
Universe and D M denotes the covariant derivative for 
a Weyl spinor. In the following, we specify more the 
definitions used here. The Dirac matrices in curved 
spacetime, T^, are defined according to the equation 
{r^r^} = — 2g fll/ . As a consequence, these matrices 
can be expressed in terms of the vierbein e° M defined 
by = e a M e^,?7afc, where rj a b is the Minkowski met- 
ric. This results in the equation = e^ a l a where 
the 7" are the Dirac matrices in flat Minkowski space- 
time. In the case of a FLRW Universe, one has <? p „ — 
a2 ( 1 i) 1 lni> which means that the veirbein are given by 
e a M = a(rj)5 a ll ,e^ a = [l/a(r])]S^ a . This implies that the 
FLRW Dirac matrices can be written as r M = [1 / a(r])]^^ . 
Since in the Weyl representation the Dirac matrices are 
anti-diagonal and can be expressed as: 



xi^ 

E M 



(25) 



where E M are the Pauli matrices in a FLRW spacetime, 
one reaches the conclusion that S M = [l/a(ry)]cr A1 where 
ct p are the Pauli matrices in Minkowski spacetime. Let 
us now turn to the expression of the covariant derivative. 
For a four-dimensional Dirac spinor 'I'd which can be 
written in the Weyl representation as 



4>a 



(26) 



the covariant derivative can be expressed through the two 
expressions: 



1 



(27) 
(28) 



The matrices a ab and a ab are defined according to the 
usual expressions, namely a ab = (l/4)(a a a b — a b a a ) and 
a ab = (l/A)(a a a b ~ a b a a ) ||. In a FRLW Universe, the 
components of the spin connection can be written as: 



W 00 = W Qia = WjOO = WijO = 0, 
w oi = w mni = 0, w j0i ■ 



(29) 
(30) 



Let us now redefine the spinor tp a by ip a = Xa/ a ( 7 ?) 3/ ' 2 • 
Then the Lagrangian C given in Eq. ( p3j ) can be re- 
written as: 
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(31) 



J a = Xad (aip) - QipdoXa- 



(36) 



This equation shows that the a spinor field in FLRW 
spacetime is conformally invariant [ p6| . Contrary to the 
case of a scalar field, there is no interaction term between 
the background geometry and the quantum field. As a 
consequence there is no creation of (massless) fermionic 
particles. Another way to put it is to say that if we start 
from the vacuum state, the system will remain in this 
state for ever. At this level, we could already conclude 
that SUSY is broken in a FLRW spacetime. Indeed the 
previous phenomenon is equivalent to say that the bosons 
possess a time dependent mass whereas the fermions still 
have a constant mass. Therefore, the condition me = 
ttif can no longer be satisfied and global SUSY must 
be broken. The phenomenon of creation of particles in 
curved spacetimes is responsible for the SUSY breaking. 

In order to demonstrate explicitly how this property 
shows up, let us pursue the calculations in more de- 
tails. Having seen that the fermions Xa behaves like free 
fermions we can canonically quantize the fermionic field. 
To do so we need to define the zero modes of the Dirac 
operator acting on Weyl fermions 



Vk + k 3 \k + k 3 



Vk Q - k 3 \k -k 3 



(32) 
(33) 



where = k\ ± ik-i- Then, the fermionic field operator 
Xa can be canonically decomposed according to: 




ik-x 



(2tt) 3 / 2 

x (V kW +(k)e lfco " + 6_, kW -(k)e- ifc °^ , (34) 

where the creation and annihilation operator satisfy the 
anticommutation relation 

{ b +M, & U<} = 5 ( k - k ')- {&-,k, bl M ,} = S(k - k'). 

(35) 

We are now in position to study the explicit breaking 
of supersymmetry due to the curvature of the FLRW 
spaces. This is the aim of the next section. 



3. SUSY breaking in curved spacetime 

The free Lagrangian in the rescaled bosonic and 
fermionic is a free supersymmetric Lagrangian apart ex- 
plicit breaking terms proportional to a' /a. The super- 
symmetric current is no longer conserved leading to a 
time dependent supersymmetric charge. A natural defi- 
nition of the current is provided by the following expres- 
sion which involves only the rescaled fields: 



Accordingly the supersymmetric charge can be expressed 
as: 



Qa = / d 3 xJ, 



(37) 



The supersymmetric charge can be expressed in terms of 
the canonical creation and annihilation operators. Using 
Eqns. (^lj) and (|34|), this leads to the relation 

dQ a _ a" f d 3 k 
dr/ a J 2ko 

+ (k)6 +ik dt k + u+(-k)& + ,_ k d 2k e lfe °" 



+u-(k)&_ >k <4 + u -(-k)6_,_ k d 2k e- jfco " ) . (38) 



Recalling that the scalar of curvature is given by R — 
6a" /a 3 we have obtained that supersymmetry is explic- 
itly broken by curvature effects. Notice that when R = 
supersymmetry is preserved. This corresponds to either 
Minkowski spacetime or the radiation dominated FLRW 
spacetime for which one has a(r)) oc 77. This results is not 
particular to the non-interacting theory but can be gener- 
alised by considering N = 1 supergravity in four dimen- 
sions. In the low energy limit when non-renormalisable 
gravitational interactions are neglected, i.e. mpi — > 00, 
the supergravity Lagrangian reduces to the curved ac- 
tion that we have considered previously. Supersymmetry 
is preserved when the background metric of the curved 
spacetime possesses enough Killing spinors, i.e. solutions 
of the spinorial equation 



D„e = 0, 



(39) 



where is the covariant derivative and e is the 
supersymmetry variation parameter. By considering 
[D^jD^Je = and the maximal symmetry of the FLRW 
metric one finds that supersymmetry is only preserved 
for flat FLRW spacetimes, i.e. R — 0. Since this is true 
even in the presence of interacting chiral superfields, this 
generalises the previous result. 

The fact that SUSY is broken in a curved spacetime 
will induce corrections to the quintessential potential 
since this one is no longer protected. We now evaluate 
the order of magnitude of these corrections. We have de- 
rived in the appendix the one-loop effective potential in 
the interacting case in the presence of curvature effects. 
The one loop effective potential reads: 



<5U eff (Q,77) 



1 



32tt 2 





1 |2 R 


2 




{ 


H --g 




ln ( 











ml 



R/6 



1 |9 R 


2 




m -1 




ln( 


_ 1 12 







m 



R/Y2\ 



, (40) 



where \i is the renormalization scale. As for the quan- 
tum corrections, the dependence in Q appears through 
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the relation m = tub = mp = d 2 V(Q)/dQ 2 . In addi- 
tion there is now an explicit time dependence due to the 
presence of the scalar of curvature in Eq. p0| ) . The evo- 
lution starts during the radiation dominated epoch where 
global SUSY is preserved since R = 0. As a consequence, 
Eq. (|4C| ) implies that 5V e g(Q,r]) = 0, i.e. no corrections 
are generated. Notice that in contrast with the quantum 
corrections the curvature effects manifest themselves at 
later times. In the matter dominated era, the effects of 
curvature no longer vanishes and the corrections modify 
the power law potential. These corrections are the usual 
quantum corrections plus additional contributions pro- 
portional to H 2 . It has been shown in the previous sec- 
tion that quantum corrections are important only deep in 
the radiation dominated era. Since H 2 is a tiny number 
during the matter dominated era, we conclude that the 
curvature effects for quintessence are therefore negligible. 



C. Kahlerian corrections 

In gauged supersymmetric models where flat directions 
are lifted non-perturbatively below a scale A where the 
dynamics of the gauge group becomes strongly coupled, 
we have seen that the effective superpotential for the am- 
plitude mode Q has the required form to lead to a in- 
verse power law quintessential scalar potential provided 
the Kahler potential is flat. The non-renormalisation the- 
orem guarantees that the superpotential does not receive 
quantum corrections due to radiative corrections. This 
is not the case of the Kahler potential which is not pro- 
tected, and therefore is modified at low energy. The low 
energy Kahler potential becomes a complicated function 
K(Q, Q*; A) which can be expanded in Taylor series as: 



K(Q,Q*;A) = QQ* 



^ Q m Q* n 

2 / a ' mn y\m+™-2 C ' C ' 

m>0,n>0 

m+n>2 



(41) 



This expansion is valid as long as Q < A which means 
that the Kahler potential will be modified at the be- 
ginning of the evolution, deep in the radiation domi- 
nated era. As already mentioned, when Q > A, the 
Kahler potential can be considered as flat. In this re- 
spect the situation is similar to what happens for the 
quantum corrections studied previously. In the case of 
the SU(N C ) gauge theory with Nf flavour there is a 
useful ansatz for the Kahler potential which illustrates 
the Kahler correction. Indeed a Kahler potential of the 
form K = (l/A f )tr(MMt + A 4 ) 1 / 2 leads to an expansion 
K w (QQt) 2 /(2A 2 ) - (QQt)4/ (8A 6) + . . . when Q <<; A 

and K « QQ^ when Q ^> A. This is typical of the situ- 
ation we are going to describe. 

Let recall the structure of the low energy Lagrangian 
in the presence of a non-flat Kahler potential. In a flat 
gravitational background the Lagrangian reads 



C = J A A x(^gdQdQ* - i 



dW 



dQ 



(42) 



where we have denoted by Q the chiral superfield whose 
scalar component is Q. The metric on the one dimen- 
sional complex curve defined by K is given by the follow- 
ing expression 



9 = 



,K. 



(43) 



Using the Taylor expansion of the Kahler potential given 
by Eq. ( pi] ) and the previous relation, one can easily 
deduce the expression of the Kahler metric: 



5 = 1 



E 

m>0,ii>0 

m+n>2 



mna mn - 



.(n-i) 



J\jn-\-n—2 



(44) 



Notice that the flat metric is modified by the quantum 
corrections in 1/A. The scalar Lagrangian can be rear- 
ranged in order to render the physical meaning of the 
non-trivial Kahler potential more transparent. Let us 
redefine the fields according to 



dQ 
dQ 



(45) 



where we are now specializing our results to the real case 
Q = Q* . The square root of g can be expanded in a series 
Jg = 1 + Xm>i &n(Q n /A n ) leading, after an inversion 
which can be carried out inductively, to: 



Q = Q 



ra>l 



A" 



(46) 



This redefinition transforms the kinetic terms of Q into 
canonically normalized ones for Q, namely d^Qd^Q. On 
the other hand, the scalar potential becomes then 



V(Q) 



1 



A 4+c 



g[Q(Q)] Q a (Q) 



(47) 



where we have explicitly shown the dependence of Q on 
Q. The metric g depends on the new field Q and is 
expandable in a series g = 1 + E n >i^K !™/A™). As a 
consequence the scalar potential becomes now an infinite 
series the expression of which is given by 



V(Q) 



A 



4+a 



E e ™ \n )' 



n>l 



(48) 



where the coefficients e n are easily computed inductively. 
This potential must be seen as the effective potential of 
quintessence in the region where Q < A. This potential 
has the usual quintessential inverse power behaviour but 
is corrected by positive powers of Q. Whether the prop- 
erties of the tracking solutions will be preserved crucially 
depends on the coefficients e n . 







Let us conclude this section by stressing that building 
a realistic model of quintessence based on global SUSY 
appears to be a difficult task. Even if these models are 
free from quantum corrections to the potential and from 
curvature corrections, the tracking properties could be 
destroyed by the Kahlerian corrections as shown by Eq. 
(fl8|). Nevertheless the main difficulty is that at the end 
of the evolution Q w mpi, rendering the SUGRA correc- 
tions unavoidable. We cannot neglect contributions in 
the scalar potential which are suppressed by the Planck 
mass. Therefore, any realistic model of quintessence must 
be based on SUGRA. The aim of the next section is to 
study such models. 

V. SUPERGRAVITY MODELS 

A. SUGRA corrections 

In this section we consider the SUGRA version of the 
model considered previously with a superpotential of the 
form W oc 1/Q a . At tree level the supergravity scalar 
potential depends on the potential G = kK + ln(K 3 |lU| 2 ) 
where K is the Kahler potential and W the superpoten- 
tial. The scalar potential is given by: 

V= ^(G i G i -3) + V u , (49) 

where the indices have been raised using the metric Gq = 
didjG and where the derivatives have been taken with 
respect to the scalar fields. The term Vb comes from 
the gauge sector of the theory and is always positive. In 
the quintessence context there is only one field and the 
Kahler potential is chosen to be flat K — QQ* . This 
leads to the following scalar potential: 

V = e«Q 2 ^(^l~(a + l) K Qi + K iQ^. (50) 

There are a few remarks at hands about this potential. 
The first term corresponds to the global supersymmetry 
scalar potential while the other two terms are supergrav- 
ity corrections in Q/mp\. There are two important effects 
of the supergravity corrections. The exponential term in 
the scalar potential introduces positive powers of Q of ar- 
bitrary degrees. Fortunately this only becomes relevant 
for Q of order of the Planck mass, i.e. for the red-shift 
z w 0. More importantly the potential can become nega- 
tive due to the second term. This implies that the model 
can become non-sensical at the end of the evolution when 
kQ « mpi. In the case a = 6, we have checked numeri- 
cally that this is indeed the case. For higher values of a 
this is still the case proving that globally supersymmetric 
models with an inverse power law superpotential do not 
resist the supergravity corrections. 

The appearance of dangerous negative contributions to 
the potential is not accidental, it stems from the —3 term 



in the potential. A possible way out is to impose that 
the scalar potential exactly vanishes and that the scalar 
potential is entirely due to a non-flat Kahler potential. 
This is what we are going to study in the next section. 

B. Moduli Quintessence 

Let us consider a supergravity model where there are 
two types of fields, the quintessence field Q and matter 
fields (X, Y % ). We assume that the gauge group of the 
model is broken along flat directions of the D terms such 
that X ^ Q,Y l = where Vd = 0. As explained in 
the previous section we impose that the scalar potential 
is positive to prevent any negative contribution to the 
energy density. This is achieved by considering that 

(W) = 0, (51) 

when evaluated along the flat direction. Moreover we 
assume that one of the gradients of the superpotential 
Wy does not vanish. With these assumptions the scalar 
potential becomes 

V = e KK K YY *\W Y \ 2 . (52) 

As expected the scalar potential is positive and be- 
comes a function of the quintessence field Q only. The 
quintessence property is achieved if this potential pos- 
sesses the run-away behaviour of the quintessence field 
to infinity. 

In the following we shall use string-inspired models 
with an anomalous U(l)x gauge symmetry [ p7|]2lj| ] . In 
the context of the heterotic string theory it is natural 
to identify the quintessence field with one of the moduli 
of the string compactification. Indeed the values of the 
moduli naturally goes to infinity. This is usually called 
the moduli problem which in fact turn into a blessing in 
the context of quintessence. Consider the compactifica- 
tion of the weakly coupled heterotic string on a Calabi- 
Yau manifold. The compactification depends explicitly 
on the moduli T a , a = 1 . . . 3, representing deformations 
of the Kahler class of the Calabi-Yau manifold. We shall 
concentrate on one moduli T a = Q on which the modu- 
lar group SL(2, Z) acts as Q — > (aQ — ib) / (icQ + d) with 
ad — be = 1 . Moreover we suppose that the gauge group 
factorises as G x U(l)x where G contains the standard 
model gauge group and U(l)x is an anomalous Abelian 
symmetry. The fields of the model split into three groups, 
the field X has a charge 1 under U(l)x and is neutral 
under G, the field Y is a matter field neutral under G 
and of charge —2 under U(l)x while the matter fields Yi 
are charged under G and possess charges qi ^= —2 under 
U(l)x- Under the modular group the matter fields trans- 
form with a modular weight — 1, i.e. Yi — > (icQ + d)~ 1 Yi, 
and similarly for Y. This is the modular weight of un- 
twisted states for orbifold compactifications. We also as- 
sume that the field X has modular weight zero. The 
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scalar potential comprises two terms: the usual scalar 
potential of supergravity and the D terms of the gauge 
symmetry. Associated to U(l)x is the D term 



V - & 
~2~ 



K X X - 2K Y Y 



QiKyYi 



ff 2 Tr(A) 

1927T 2 ' 



( P1 



(53) 



where gx is the U(l)x gauge coupling and g the unified 
string coupling. Modular invariance is compatible with 
the Kahler potential 



K = XX* - -In [Q + Q* - k\Y\ 

K 



(54) 



The D term potential vanishes altogether along a flat di- 
rection where the field X acquires a vacuum expectation 
value (vev) breaking the Abelian symmetry U(l)x at 



(X) = 



-ff 2 Tr(X ) 

1927T 2 



-mpi, 



(55) 



while the other fields vanish altogether. The value of the 
vev (X) is equal to a scale f which is fixed in the heterotic 
string. This is not the case any more in the context of 
type I strings where £ is a moduli which is not fixed in 
perturbation theory. Expanding the superpotential in 
terms of Yukawa couplings, we obtain 



W = \{Q)X 2 Y + 



(56) 



where we have only taken into account the coupling such 
that Wy = X(Q)(X) 2 along the flat direction. The func- 
tion X(Q) is a modular form of weight —2. The scalar 
potential along the flat direction is 



V(Q) oc 



W 4 A 2 (1) 
12Q 2 



(57) 



As already mentioned for the Kahler corrections, the 
kinetic term in such models is not standard. In the 
present case, the kinetic term for the moduli Q is 
[3/(4k<5 2 )](9Q) 2 . Therefore, it is more convenient to 
redefine the fields such that the kinetic term becomes 
standard. This is achieved by means of the following 
expression: dQ — \/3/{y2KQ)dQ. The scalar potential 
then transforms into: 



V(Q) 



A 2 (l) 
12 



(X) 



(58) 



It has been shown in Refs. [ |1 8|j2 S| ] that this potential also 
possesses remarkable tracking properties even if it has no 
inverse power factor of the quintessence field. However it 
suffers from phenomenological problems. First of all the 
slow decrease of the potential and the large value of the 
vev (X) imply that the quintessence field will have to be 



much larger than the Planck scale at the end of its evo- 
lution. This is a drawback as this would be in the string 
regime where the supergravity approximation is not valid 
anymore. In addition, one can show that the equation of 
state for this potential is such loq — lub both in the radi- 
ation and matter dominated epochs J2j| . The scalar field 
follows exactly the behaviour of the background. As a 
consequence, the value of £Iq has been shown Jl^,^] to 
be limited to the relatively low value fig w 0.15. This is 
not enough to reproduce the data which seem rather to 
indicate that Oq ss 0.7. Nevertheless, the fact that the 
moduli field possesses a tracking behaviour is relevant as 
it could represent one of the components of the energy 
density of the universe. 



C. An inverse power law SUGRA model 

We now present a toy supergravity model with an in- 
verse power law potential. To our opinion, this model is 
the most interesting one although we only present it as 
an existence proof of the quintessence property in super- 
gravity. 

Let us use the same framework as in the previous sec- 
tion. We assume that the U(l)x Abelian symmetry is 
broken by the Fayet-Iliopoulos D term. The superpoten- 
tial is expanded as 



W = XX Y 



(59) 



where A is taken to be constant. This superpotential 
preserves the gauge symmetry of the model. As claimed 
earlier we do find that (W) = along the D-term flat di- 
rections. The main difference between the present model 
and the case of the moduli field of the heterotic string 
is the different form of the Kahler potential. We choose 
the geometry of the moduli space to be singular at the 
origin, namely 



K = XX* 



(QQ* 



2p-2 

1 1 



\Y\ 



(qq*y 



(60) 



where m c ~ 10 15 GeV is the cut-off of the theory. It 
is reasonnable to choose it of the order of the unifica- 
tion scale. Higher order terms in \Y\ 2 can also be in- 
cluded. This is the only relevant terms in the Kahler 
potential if one assumes the existence of a modular sym- 
metry. This modular symmetry is also important to pre- 
vent any Kahler correction. The scalar potential of this 
supergravity model is given by 



V(Q) 



A 4+a 
Q a ' 



(61) 



where it is understood that Q is now the canonically 
normalised field. The constant A is given by A 4+Q = 
2 n / p X 2 (X) 4 m" where a = 2n/p. This potential has been 
studied in details in Ref. J]j|. There, it has been shown 
that, despite the appearence of positive powers of the 
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field, the tracking properties are completely preserved. 
Using the constraint (X) > 10 2 GeV as the Abelian 
symetry has to be broken above the weak scale we find 
that the fine tuning problem can be overcome provided 
that a > 11, see Ref. |lj|. This means that in order to 
obtain J7q w 0.7 today, one should fix A to a value which 
is not far from the natural scales of high energy physics. 
For a = 11, actually we have A ss 10 11 GeV. The evo- 
lution of the energy density of the quintessence field is 
given in the following figure. 




5 10 15 20 25 
log (1+z) 



FIG. 5. Energy density in the SUGRA model of 
quintessence for a = 11. The solid line is the energy den- 
sity of quintessence whereas dashed-dotted and the dashed 
lines are the energy density of radiation and matter respec- 
tively. The dotted line is the energy density of quintessence 
for the potential V{Q) = A 4+a Q' a for a = 11. The initial 
conditions are such that equipartition is realized just after 
inflation. 

The evolution of the equation of state loq for this model 
is displayed in the following figure. The evolution of the 
equation of state is almost unchanged during all the cos- 
mic evolution. This is due to the fact that during this 
epoch Q/mpi is very small. As a consequence the ex- 
ponential factor in the SUGRA potential plays no role 
and the SUGRA potential reduces to the usual poten- 
tial. The situation changes at the end of the evolution. 
Since the field is on tracks one has Q w mp\. This time 
the exponential factor in the potential plays a vital role 
and modifies the value of loq today. 
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FIG. 6. The dotted line represents the evolution of loq for 
the potential V(Q) = A 4+a Q~ a with a = 11 whereas the 
dashed line represents the evolution of the equation of state 
for the same value of a and for the SUGRA model presented 
in this article. 

We see that the value of the equation of state today, 
for a = 11, is given by Jl5| : 

loq w -0.82. (62) 

This is a remarkable property since in the context of usual 
tracking solutions it is not possible to obtain a number 
less than ss —0.7 |ll[] whereas the observations seem to 
indicate that the value of the equation of state today is 
rather such that — 1 < loq < —0.8. We would like to note 
that this nice property is almost independent of a. This 
is illustrated on the following figure where the relation 
loq — a is displayed 

0.0 L 

-0.2 ; : 

-0.4 '- 1 

at 

3 

-0.6 ; ~ 
-0.8 

-1.0 \ ................. ■ 

5 10 15 20 

a 



FIG. 7. Equation of state loq in function of a for the 
SUGRA model of quintessence proposed in this article. 

The equation of state is almost independent of a be- 
cause its value today is roughly speaking determined by 
the exponential term in the potential which is a indepen- 
dent. As a consequence, no fine tuning of a is required 
in order to obtain a reasonable value for loq. 



12 



The implications for structure formations of this model 
seem to be also very interesting and are currently under 
investigations |3(J. 

More details on this model can be found in Ref. Ila] . 



VI. CONCLUSIONS 

We have studied the quintessence scenario recently 
proposed in Ref. (r|. We have concentrated on mod- 
els with an inverse power law potential and the diverse 
corrections induced by its embedding into High Energy 
Physics models. More particularly we have focused on 
the quantum corrections in the non-supersymmetric set- 
ting, the curvature and Kahler corrections in the super- 
symmetric case, and the effect of the non-renormalisable 
interaction suppressed by the Planck mass in the super- 
gravity context. We have verified that quintessence is sta- 
ble to the one-loop quantum corrections, preserving the 
existence of tracking solutions at small red-shift. We nev- 
ertheless argue that the solution to the fine-tuning prob- 
lem requires to consider the quintessence models within 
the framework of particle physics beyond the standard 
model. Emphasizing the supersymmetric scenario, the 
most likely candidate to describe the physics beyond the 
weak scale, we generalize the usual inverse power law 
potential of the SU(N C ) super-QCD case to more gen- 
eral supersymmetric gauge theories with Dynkin indices 
fJ- < Madj ■ In supersymmetric models the quantum correc- 
tions to the superpotential vanish. However there are two 
types of potentially dangerous corrections. Firstly the 
smallness of the effective mass of the quintessence field 
implies that one must consider the small SUSY breaking 
induced by curvature effects. We show that these effects 
are not a threat to the quintessence scenario. Secondly 
we investigate the corrections to the Kahler potential and 
show that they play only a role at the beginning of the 
evolution, decoupling at a scale of the order of the gaug- 
ino condensation scale, and are strong enough to jeopar- 
dize the quintessence property. Far more relevant is the 
necessary inclusion of dangerous supergravity corrections 
due to the presence of Planck mass suppressed interac- 
tions which induce negative contributions to the potential 
at small red-shift. They destroy the quintessence prop- 
erty if the quintessence models are not treated within the 
realm of SUGRA. 

This is what we do in the final section considering 
models where the expectation value of the superpoten- 
tial vanishes altogether. This prevents the appearance of 
the negative contributions to the energy density. We give 
two explicit models where this scenario is at work. One is 
based on the heterotic string at weak coupling, the role of 
the quintessence field being played by a moduli. The run- 
away behaviour of the moduli is of the exponetial type 
limiting the possible contribution to the energy density 
to Qq = 0.15 per moduli. In usual orbifold compacti- 
fication this can amount to a Qq = 0.45 contribution. 



Finally we present an existence proof of an inverse power 
law potential in SUGRA by constructing a model based 
on a particular Kahler potential. This model is partic- 
ularly promising as it provides a scenario for which the 
value of ojq is —0.82 for Q m = 0.3. This is within one 
sigma of the recent experimental analyses Q . 

Let us emphasize that the models of quintessence 
within the framework of SUSY or SUGRA suffer from 
a problem raised in |T^ | concerning the necessary break- 
ing of supersymmetry. Indeed there arc two aspects to 
this problem, both linked to the types of corrections con- 
sidered here. Let us first deal with the SUSY case. As 
advocated in p2j the present of a non-zero F = M| term 
of the order 10™GeV leads to an intolerably large value 
of the cosmological constant. This necessitates to invoke 
the prejudice already used in our paper that one does not 
know the mechanism which would force the cosmological 
constant to vanish. Quintessence scenarios do not aim 
at solving the quantum cosmological constant problem, 
see the conclusion of Ref. jl8| . Another difficulty springs 
from the possible non-flatness of the Kahler potential and 
a coupling between the quintessence field Q and the field 
responsible for the supersymmetry breaking leading to a 
term given by: 

SV(Q) = k(Q)Ml, (63) 

where the series expansion of A: is a power series sup- 
pressed by the Planck scale. It gives a large contribution 
to the mass of Q and the cosmological constant. 

Of course one should repeat the arguments within the 
SUGRA context. At tree level one can always fine-tune 
the minimum of the potential to be zero. The power law 
expansion of the potential might just give a correction to 
the quintessence behaviour which would only modify the 
end of the evolution of Q while preserving quintessence. 
In the explicit model of section V.C, the SUSY breaking 
is induced by the F terms of the dilaton and the moduli. 
The dilaton does not couple directly to the quintessence 
field Q in the Kahler potential because of the anomalous 
[/(l)x symmetry. This implies that our model is not 
affected by a dilatonic SUSY breaking. Assuming the 
existence of a modular symmetry, we find that the effect 
of the SUSY breaking by the moduli is to induce terms 
like F 2 /Q 2p in the scalar potential. More studies are 
required to evaluate the effect of this term, in particular 
its order of magnitude. However, the fact that this term 
is proportional to the inverse of the quintessence field 
suggests that it will cause no problems. This question 
will be adressed elsewhere [ f?l| . 

In conclusion we would like to emphasize the explicit 
construction of particle physics models leading to the 
quintessence property is still to be further developped. 
It seems to us clear that the indications given in this 
paper show that the building of SUGRA models should 
turn out to be very fruitful. 
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VII. APPENDIX I: THE EFFECTIVE ACTION IN 
CURVED SPACE 

In the following we shall consider the case of curved 
spacetimes R ^ 0. Let us consider the following La- 
grangian 



m(w) , , 
+ — 7T-W> + 



fh(<p) 



tpi) + V(<p) 



(64) 



where, for convenience, we have not written the spinorial 
indices. The mass and the potential are related to the 
superpotential W as 



i(if) 



d 2 W 
dip 2 



V(<p) 



dW 



dtp 



(65) 



This is the usual globally supersymmetric Lagrangian 
coupled to gravity. As supersymmetry is not preserved 
by the background curved space the Lagrangian receives 
quantum corrections. The effective potential is renor- 
malised due to the background geometry. Let us calcu- 
late this potential at the one loop level. 

It is technically necessary to perform a Wick rotation 
to obtain a Euclidean field theory on a curved Rieman- 
nian manifold. This is globally feasible if one chooses 
the curved spacetime manifold to be a globally hyper- 
bolic manifold with a Lorentzian signature. This implies 
that the time coordinate is globally defined allowing us 
to perform the appropriate rotation t — > it. In order 
to have a good control on the Feynman path integral 
it is convenient to expand the bosonic field ip around 
the classical solution ip c of the Klein-Gordon equation 
— Aip c + (dV/dip*) = where A is the curved Laplacian 
on the Riemannian manifold. Let us denote ip — ip c + dp 
where dp is a full-fledged quantum field. The effective po- 
tential does not involve derivatives of ip c , we will only 
retain non-derivative terms in the expansion of the effec- 
tive action 



-Seff(Vc) 



VipVipVipe- 



where S c g is defined by the equation 
Seff = J <^X\fgV cii {ip c ). 



(66) 



(67) 



The effective potential is the sum of the classical part 
corrected by a quantum contribution 



V eS (ip c ) = V(ip c ) + 6V(ip c ). 



(68) 



The correction term is due to the integration over the 
Weyl fermions and the bosonic fields. To leading order 
the effective action is 

Ses&c) = trln(A + \m\ 2 ) - Jtrln(A F ), (69) 



where Ap is the operator Ap = D F Dp and 
, iD m 



(70) 



where D and D are the Dirac operators acting on Weyl 
fermions of both chiralities. The previous expression can 
be computed using the £ function regularisation where 
the determinant of an operator is det A = exp[— Ca(0)] 
and the ( function is ((s) = 5Z n l/A* as a function of 
the eigenvalues A„. The (a function is related to the 
heat kernel solution of 



d 

AG(x, x', a) = ——G(x,x',a), 



(71) 



with the boundary condition G(x, x',0) = S(x — x') via 
the Mellin transform 

Ca(s) = J™ / d A x^G(x, x, a). (72) 

The scale /x is a renormalisation scale. The determinant 
of the operator A is essentially given by the asymptotic 
expansion of the heat kernel 



G(x, x, a) 
leading to 



1 



(47TC7) 2 



1 -Ticr + 0(cr 2 ) 



(73) 



det^ = ex P (/d^^(m^-^)). (74) 

The asymptotic expansion of the Heat kernel of the op- 
erator —A + \m\ 2 gives t\ = —R/6 + \m\ 2 and for the 
Dirac operator it yields rf — \m\ 2 + R/12 where R is the 
scalar curvature. The one loop effective potential is then 



SV eS (Q,r]) 



32tt 2 





1 |2 ^ 


2 




{ 


H - g 




ln^ 











R/6 





2 








H 









\m\ 2 + R/12 
J 2 



(75) 



Notice that the one loop correction vanishes if the curva- 
ture is zero. 



VIII. APPENDIX II: SUPERGRAVITY VACUA 

In this appendix we shall be interested in finding 
whether a given FLRW space time preserves supersym- 
metry. The FLRW background does not break supersym- 
metry in an explicit manner when the gravitino ip^ = 
is invariant under supersymmetry transformations. We 
have assumed that the quintessence field Q is only cou- 
pled to other field via the gravitational interactions. In 
a hidden sector the supergravity can be broken with a 
non- vanishing mass m 3 / 2 leading to the condition 
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D^e 



1-3/2°^ 



0, 



(76) 



where e is a Weyl spinor representing the supersymme- 
try variation. Computing the commutator [D M , D„] and 
using R^pa = (fl/2)(r? w r? !/tr - rj^pTj^), we find that 
the supersymmetry is generically broken apart from two 
cases. If the gravitino mass vanishes, the supersymme- 
try is fully preserved by FLRW spaces with vanishing 
curvature. When the gravitino mass is not zero the su- 
persymmetry is not broken if the FLRW space is the 



anti-Dcsittcr manifold with curvature R = 
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